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?ITLE 



TYPE 

NUMBER OF WORDS 

TEMPORARY STORAGE 

DURATION 
ACCURACY 



PARAMETERS 



1.7 Precision Floating Binary Arithmetic and Double 
Precision Fixed Point Arithmetic With Floating (DOI or SADOl) 
Decimal Conversion 

Interpretive routine entered like a closed subroutine a 
610 

Locations 0, 1, 2 and 6 locations specified by a preset- 
parameter S3, S3-5S5- 
See the order codec 

Floating point numbers are rounded to 68 binary places 
(together with a sign) . Fixed point numbers are carried 
to a full 78 binary places (together with a sign)o Print 
out is rouadedo 

S3| during input of the program location 3 must contain 
c -39 



t x 2 " where t is the location of the first word of 

temporary storage 
DESCRIPTION This routine was written as a flexible general purpose 

double precision routine e It Is suitable for problems requiring twelve tc twenty- 
three decimal places of accuracy which do not require a large amount; of computing 
time* As an example of problem sizes inversion of a 30x30 matrix takes approximately 
kO minutes using an interpretive matrix inversion routine 

The method of Interpretation and the interpretive, order 
code are taken almost directly from ILLIAO library routine Ale It is suggested 
that a potential user of this routine who has not used either Al or this routine 
should also read the description for routine Al, Coding methods and examples are 
discussed there „ 

Because this routine has some special features it is longer 
than a similar routine, ILLIAC library routine Ak - It is slower or approximately 
the same speed as routine Ah- in its arithmetic operation «, The special features 
referred to are 
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(1) a closed subroutine for standardizing numbers which is very- 
fast 

(2) closed subroutines for multiplication, division, sign change, 
read a fixed point fraction or integer from tape, print an 
integer = These subroutines can be used by auxiliary sub- 
routines * 

(3) Fast input which reads numbers at maximum reader speed 
(k) Ability to read numbers punched in the form sign, integer 

part, decimal pointy fractional part. 

(5) Dual fixed point - floating point operations with changes 
from one mode to the other mode by execution of interpretive 
instructions . 

(6) Interpretive instructions for adding and subtracting B-registers * 

(7) ^096 addressable numbers (pairs of memory locations some of which 
are on the drum) a 

(8) In floating point division the maximum error is less than *52 in 
the least significant quotient "bit- If the divisor should divide 
exactly, the quotient is exact. 

NUMBER STORAGE, METHODS OF NUMBER REPRESENTATION 

Locations S3, 133* 2S3 will be called the floating 
accumulator j locations 3S>, 4S3> and 5S3 will be called the number register . 
In preparation for each arithmetic operation an operand is placed into the 
number register from the memory. The arithmetic result is placed in the floating 
accumulator. Numbers are not brought from the memory before the execution of 
non -arithmetic instructions , 

During floating point operation a number is stored in the 
memory in a packed form requiring two consecutive memory locations n, n+1. It 
has the following standard form 

n 4 — -X-,— — ^ n+1 
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-39 E x 
and x = (x, + 2 x p )*2 .In the standard form either 

x, = x = and E = -512 
12 x 

-39 / 

or one of -1 < x + 2 ^ x < - 1/2 

1/2 < x n + 2~ 59 x < 1 

with -512 < E < 512 , Thus, of 80 digital positions , the exponent uses 10, 

X / -39 n 

the positive sign of the least significant part L(x 1 +2 x ) uses 1, and the 

-39 r 
fractional part x + 2 ' x uses 69 digital position. 

During fixed point operation numbers are stored in an 

unpacked form in locations n, n+1, in the following way* 




n+1 I x 



When a number is brought from the memory and is put into 
the number register, it has the . following form? 

N(5S3) = x x 
N(^S3) = x 2 

N(5S3) = 2~ 55 (E +512) floating point 
= fixed point 

When a floating point number is unpacked the ten digital 
positions of ^S3 which hold the exponent are cleared to zero. 

Arithmetic floating point results in the floating 
accumulator have the form 

2[B(33) + 2" 39 N(1S5)] - 2^ H(2S3)-512] = a „^ 

Results in the floating accumulator are not standardized after arithmetic . 
Standardization may be accomplished by use of the interpretive N2 F instruction 
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Arithmetic fixed point results in the floating 
accumulator have the form 



2[N(S3) + 2" 39 N(1S3)] 2 2 W 2S 3)J = 



answer 



In this case, N(2S3) > and N(2Sj) > results from the automatic right 
shift of overflowed numbers resulting from addition, subtraction and 
multiplication. To perform a fixed point division using this routine the 
programmer must guarantee division overflow will not occur. The division 
subroutine further requires the divisor y to satisfy 

|y| > 1/2 . 

NUMBER STORAGE, ADDRESSING OF LOCATIONS 

The 1024 locations in the Williams Memory together with 
additional locations on the drum are directly addressable using a 12 digit 
address carried by an interpretive instruction together with use, if desired, 
of 12 digit b-modifierse A 12 digit address plus a 12 digit modifier can be 
used to construct any address in the range < address < 8l91„ The 
correspondence between interpretive/ addresses and ILLIAC addresses are 

Interpretive - 1023 Williams Memory - 1023 
Interpretive 1024 - 8191 Drum 2560 - 972 7 

As an example, if the interpretive address (unmodified) 
is 1023, then the most significant part is at Williams Memory location 1023 • The 
least significant part is at drum location 2560. 

INTERPRETIVE INSTRUCTIONS 

The interpretive . instructions have the same form as the 
regular ILLIAC instructions. This consists of two function digits (two sexadecimal 
characters) followed by a single 12 -bit address ■. 

The first sexadecimal function digit designates a b register 
to be used in executing the instruction. The correspondence is as follows; 
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A first function digit of 0, 1, . .«, 7 refers to b-registers 
0, 1, •«», 7 respectively- 

A first function digit of 9, K, S, «•»<, L refers to "b-registers 

1* 2, Js * • • j» ( * 

A first function digit of 8 does not refer to a b-register. 

The second sexadecimal function digit indicates the 
type of operation to he performed. This is discussed in detail below. 

The actual structure of the b-registers will be 
described in the sequel. However, each b-register contains two parts, an 
addr es s -register and a count register. The former contains 12 binary stages 
while the latter contains 20 binary stages = 

Let F denote the number in the floating accumulator and 
F(n) the number in locations n, n+1* Let 11^,11,, .-.., n„ denote the addresses and 
c s ■ c 9 »..,- c 7 the counts in b registers 0, 1, »*.* 7 respectively. In 

the sequel numbers F(n+n ) will be indicated using the following conventions! 

13 
no = and n + n is taken mod. 2 . 

Interpretive instructions may be stored only in the Williams Memory * 

bO n < b < 80 Form F - F(n + n,) 

bl n .0 < b < 8 Replace F by - F(n + n. ) 

b2 n < b < 7 Replace c, by c + 1, n, by el + 2 and jump to the 
interpretive instruction at the right side of location n 
(0 < n < 1023) if c + 1 > 0. If c + 1 < take the next 
interpretive instruction «, 

82 n Test the sign of the floating accumulator '. If F > jump to 
the interpretive instruction at the right side of location n 
(0 < n < 1023)<. If F < take the next interpretive instruction* 

92 n Jump to the interpretive instruction at the right side of location 
n. (0 < n < 1023) 

K2 Change the mode of operation from floating point to fixed point. 

S2 Change the mode of operation from fixed point to floating point. 

W2 Standardize the floating accumulator » 



Instructions 92, K2, S2, N2, 93, K3, S3, N3 do not effect any change in the 
b -registers but they do make a reference. These variants might be used before 
the 8L instruction ■ 

b3 n < b < 9 • 

Same as corresponding b2 n instruction except that the jump 

is to the left side of the location. 
b3 b = K, Sy N. 

Same as corresponding b2 instructions . 
b4 n < b < 8. Form F + F'(n + n) 
b5 n < h < 8. Replace F by F(n + n) 
b6 n < b < 8. Form F / F'(n + n. ) 
b7 n < b < 8. Form F x F(n + n) 

88 Eead one number from the input tape and put it in the floating 

accumulator. See the description below of the punching of data 
tapes . Numbers on input are converted to floating binary numbers 
and will not be standardized unless they are stored. Floating 
decimal numbers may have 2k or less decimal digits in the fractional 
part . 

89 n 1 < n < 22. Punch or print F as a sign, one space, n decimal digits,- 

one space, sign of exponent, and the exponent to 3 decimal digits 
with zero suppression on the first two digits . F is converted before 
punching or printing to a floating decimal number. The floating 
accumulator is not standardized. This instruction destroys the 
contents of the floating accumulator. 

Format may be controlled by use of the 8f n and S>¥ n 
instructions which control the punching of CR-LF characters and 
spaces between groups of decimal digits* 

The bO, bl, b4, b5, b6^ b7 instructions with 
b = 9,> K, S, . . ., L do the same thing as the corresponding instructions with 
b = 1, 2, 3, ..., 7- 

The b8 and b9 instructions with b -f- 8 do the same thing 
as the 88, 89 instructions but in addition make a reference to a b-register. 
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bK n < b < 7. Set n = and c = -n. 
bK n b = Q, 9, ... y L. 

Used to place an integer in F. If n < 512 set F = n in 
standardized form. If n > 512 set F = n-1024 (standardized), 
bS n < b < 8 Replace F'(iW-n ) by F£ 

The number F is standardized before being stored at location 
n +-n, . The contents of the floating accumulator are not 
changed * 

When the binary exponent E of F becomes too large or too 
small the following special conventions apply. E < -512. The number stored at 
location n + n, has fractional part = 0, exponent = -512- The machine representation 
of E is E + 512 so that underflow results in two ILLIAC zeros being stored. 

E > 512.o ILLIAC stops on FF 039 • This stop occurs only during the execution of the 
bS n instructions 

bN n < b < 8 Replace F by |f | - |F(n + n, | 

8 J n Leave the interpretive mode of operation and jump to the ordinary 
ILLIAC instruction on the left side of location n. To re-enter the 
interpretive mode of operation after an 8j n instruction one may 

(1) use a normal subroutine entry to A7 

(2) jump to the left instruction at 17 L of A7» In this case 
the next interpretive instruction is the instruction 
following the 8j n instruction. 

(3) use other special purpose entries described below. 

bJ n < b < 7, Let m be the integer whose digits are the least 

significant 12 binary digits of the number at Williams Memory 
location n* Replace il by a-m (mod 2 ^) and leave c 
unchanged. This instruction is especially useful when n is 
the address of a b-register because this instruction then 
differences 2 b-registerso 

8F n If n > 0, punch or print one CR-LF nowj one CR-LF and one 2 -hole 
delay before the next 89 instruction^ an CR-LF and delay before 
every nth such 89 instruction thereafter. For example if n = 2 a 
CR-LF and delay are punched before the first, third, fifth, ... 89- n 
instruction executed. 



If n = Oj one CR-LF character is punched or printed "but 
otherwise the instruction acts as if n = k-096o 

After execution of each 89 n instruction a 5 hole 

delay and two spaces are punched or printed unless a CR-LF will precede the 

next 89 n instruction 

9F n During the punching or printing of all succeeding numbers } 

punch or print one space after each n digits of the fractional 
parte 

If n = 0, or y if this instruction is never executed , 

no spacing occurs . 

12 
bF n ■< b < 7* Same as bj instruction except that ri + m (mod 2 ) 

is formed . Does not change c, • 

• 12 
bL n < b < 7° Replace ri by il + n (mod 2 ) Does not change 

v 

bL n b = 9j K, ... ; X. Replace n, } n , . .., n respectively by 

12^ 
n., -n. n -n. . ., n„-n (mod 2 ). Does not change c. . 
1*2 7 b 

8L n If b is the number of the last b-register referred to, replace 
n, by n. Does not change e_ . 

OPERATION TIMES The following table of approximate operating times is 
offered as a guide to the programmer . The time for each instruction involves 
three component times,, interpretation, memory access, and time of execution . 
Further, the drum acts as a slow access memory for which only a random access 
time can be given . 

The time for interpretation of the digital combinations 
in an instruction is so nearly the same for b = 8 vs b £ 8 that no distinction 
is made between the two types of instructions . 

Interpretation; 1.6 ms per instruction 

Memory access before bO, bl, bk 9 \fj 9 b6, b7, b.N 

instructions % TT -m- *« o « 

Williams Memory 2 °2 ms 

Drum lk-°k ms random access 



Times of execution 

bO 3°k ms 

bl 1.4 ms 

b2 r 3 < b < 7 1-5 ms 

82,3 1.0 ms 

92,3 1.0 ms 

h4 2.5 ms (generally) 

1.0 ms for F + 

1.8 ms for + F(n + a ) 

b5 .9 ms 

b6 7.9 ms 

b7 4.9 ms 

88 3-5 his per character + 30 ms per number 

89 n Minimum time 176 + l6n ms* 

If 9F instruction is used, add l6 ms for each space 
■which is punched. For example, a 20 digit number with 
spaces after every 3 digits, 

176 + 16 x 20 + 16 x 6 = 592 ms. 

bK < b < 7 .7 ms 
8k n n=0 .7 ms 
8K n n/0 3-7 ms 
bS Williams Memory 5.2 + .Ol85(3k + 2k [k/8]) ms where k 

is the number of shifts to standardize and [k/8] the 

greatest integer < k/8. 

Drum Memory 17.O + .Ol85(3k + 2k [k/8]) ms random access 
bj < b < 7 1.3 ms 
8j .2 ms 

bF < b < 7 1.3 ms 
bL < b < 7 1.0 ms 
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EXAMPLE OF TIME ESTIMATION 

The following is a loop which can he used to compute 
the inner product of a sequence of numbers stored in the Williams Memory with 
a sequence of numbers stored on the drum 

. OK nF 
8K F 

1 . 8s kF 

05 2000F 

2 . 07 900F 

Qk ky 
3° 8s ky 

02 1L 

In this loop 5 instructions are interpreted, requiring '5 x 1.6 = 8.0 ms . Two 
Williams Memory accesses require k.k ms. Execution times are 

o9 + ko9 + 2*5 + 5.2 + 1.5 = 15.0 

Total time exclusive of the drum access, 

8,0 + kok + 15.0 = 27A ms. 

Consequently the loop time is the time of 2 drum revolutions and 5 sector times 
or about 35 ms * 

ENTRIES TO A7 The closed subroutine entry to A7 is 

p 50 pF 

p + 1 26 qF 

where q is the location of the first word of A7 = The first interpretive 
instruction must be contained in the right half of location p+1. 

After leaving the interpretive mode by use of the 8j 
instruction, A7 may be re-entered in several ways . A closed subroutine entry 
may be used. A jump to the left side of 17L (relative to the first word of A7) 
results in the interpretive instruction following the 8j instruction being 
interpreted next, A jump to the left side of 2L will cause the 8j instruction 
to be executed again. A jump to the left side of 90 L causes standardization of 
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the floating accumulator, then preceding to the next interpretive instruction . 
This assumes the floating accumulator is not overflowed. A jump to the right 
side of 139 L causes a test and correction of the floating accumulator for 
overflow. The next interpretive instruction is then taken. To use this entry 
Q must agree with 1S3 • A jump to the left side of l^k L clears the floating 
accumulator. The next interpretive instruction is taken » 

B-KEGISTERS The b -registers are stored in eight consecutive 

locations 38L, 39L> °-, ^L of A7- If we denote a word by a Q , a , . .. a then 
a„ ... a ir . constitute the count c, carried by the b-register, a or ,.».a__ are 

U x.y D 2.U d. ( 

always zero° a o«* ° * a ^o constitute the address tl carried by the b register. 

The programmer should remember that storage of a number 
requires two locations Consequently the b2 n and b3 n instructions, the 
count-jump instructions, form il + 2 . 

Instructions of the form b,F 3° + b L will cause the 
addition il + zl = il* . Instructions of the form b J 38 + b L will cause the 

subtraction il - n = n* . 

PUNCHING DATA TAPES FOR INPUT 

Numbers may be punched on data tapes as floating decimal 
numbers or as a sign, integer part, fractional parte In either case, the number 
when read is converted to a floating binary number which is put in the floating 
accumulator. It will not in general be standardized. No special instructions or 
markings are necessary to distinguish the two cases. In the fixed point mode, 
numbers with zero exponents are input to full double precision accuracy » The routine 
determines the punching format from the 5-hole characters read; each number on 
tape must be terminated by a non-space character « 

Floating decimal numbers having 12 or less decimal digits 
in the fractional part may be punched as follows 

1. Sign of number 

2* decimal digits of the absolute value of the fractional part 

3. sign of the exponent 
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ko decimal digits of the absolute value of the exponent. 

Non significant zeros on the left may be omitted. 

For example 052 or 52 » 
5 . Any 5-hole character other than a space causes reading 

of exponent digits to end . 

The input routine is sensitive to all 5 hole characters except spaces. Spaces 
may be used for punching out errors . So long as the number of decimal digits 
in the fractional part is 12 or less, one 5-hole character other than a space 
or decimal point may be used anywhere between the 1st decimal digit of the 
fractional part and the sign of the exponent. Consequently output tapes by 
library routine Al may be read by routine AJ } but not conversely . 

Floating decimal numbers having 13 to 2k decimal digits 
in the fractional part must be punched as 

1. Sign of the fractional part* 

2. 1s~ digits of the fractional part, k < 12, a 5 hole character 
other than a space or decimal point, k p digits of the 
fractional part, k p < 12. 

3-5= As for numbers with less than 13 digits in the fractional 
part . 

Spaces may be used, again, for format or error correction. When 12 or more 
decimal digits occur in the fractional part, there must not be any other 5-hole 
characters used except as described above. 

Numbers in the integer-fraction representation are 
punched as 

1. sign of the number 

2 . integer part of the absolute value 
3« decimal point 

k. fractional part of the absolute value 
5* any 5-hole character other than a space. 
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If either the integer or fractional part is zero it need not be punched . Non- 
significant zeros may be omitted to the left of the decimal point, trailing 
non -significant zeros to the right of the decimal point may be omitted. Hence 

= +o =-.=+.0 

1 = +1. 
l/lO = + .1 • 

The integer part can have 12 or less decimal digits, the fractional part 12 or 
less decimal digits <> 

The integer-fraction representation was included as a 
convenience for the manual preparation of data tapes . The routine has no 
corresponding output facility. Further, where the greatest possible accuracy 
is desired for the floating binary number put into the floating accumulator, the 
floating decimal representation should be used on data tapes . During reading of 
a integer -fraction, the integer part is read and stored in S3, the fractional 
part is read and stored at 1S3* The exponent 39 is assigned . If for example the 
number is +.1, this is clearly less accurate than computing l/lO by a double 
precision division, the process used when reading floating decimal numbers . 

The output of the routine can, of course, be read by the 
routine, and hence conforms to the above rules . A punching error on a tape 
which causes the above rules on format to be violated results in an FF 03K stop. 

The following descriptions are primarily intended for 
those who wish to write fast auxiliary subroutines <, 

STANDARDIZATION SUBROUTINE 

The subroutine is entered using a pair of instructions 

P 1 F5 pF 
p + 1. I 26 2&%L 
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"where the address of the jump is relative to the 1st word of A7° The most 
significant part of the double precision number should be in 3S3> and the least 
significant part in Qo When the routine is left,, the standardized number is 
in AQ. The number of shifts used in standardizing the number is in 5S3* If the 
number was zero, 5S3 is set to zero. Except for zero the standardized number 
satisfies the inequalities l/2 <x<lor-l<x< -l/2 . Temporary storage, 
3S3, 4S3, 5S3- Operating time 1.6 + <Ol85 (3k + 24[k/8] ) ms, [k/8] denotes the 
"greatest integer less than"o 

MULTIPLICATION SUBROUTINE 

This subroutine multiplies the double precision fixed 
point number in S3., IS 3 by the double precision fixed point number at n, n+1. 
It is entered by the instructions 

p j\j hi- j. jr s j px 1 



p + 1 26 230L 

When the subroutine is left, the rounded double precision product is in S3., 1S3' 

The least significant part is also in Q= Temporary storage, location 1. 

-39 -39 

If the two numbers are x +2 x p and y + 2 "^ y 

the "product" which is formed is 

--39 /„ „ . „ „ , 2/ll „ n -39- 



x^ + 2 jy (x x y 2 + x 2 y 1 + 3 A x 2~ jy ) 



This product is truncated to a double precision number » The rounding factor 
3/k x 2~^" is justified on the basis of l/2 + expected (x y ) = l/2 + l/k = 3/h 

DIVISION SUBROUTINE 

This routine divides the double precision fixed point 
number in S3* 1S3 by the double precision fixed point number in 3S3* Q° It is 
entered by the instructions 

p F5 pF 



p + 1 26 244L 
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When the routine is left, the double precision quotient is in S3> 1S3. Temporary 
storage locations used are 0, 1, 3S3, ^-S3- 

The method of division is as follows i The sign of the 
divisor Is stored at location 0. The absolute value of the divisor is computed 

and stored in 3S3* ^S3° 

-39 -39 

Let x + 2 x p be- the dividend and y + 2 y y be the 

absolute value of the divisor. First a quotient q.. satisfying 

x x + 2 Jy x 2 = q x y x + 2 Jy r, y 1 > r > 
is computed. The quotient q and remainder r are then corrected so that 

x x + 2~ 59 x 2 = q£(y x + 2~ 59 y 2 ) + 2"" 59 (r* + (^ - q*) y 2 ) 

with y, > r* > 0. This correction results in a double precision r* A correctly 
rounded quotient r*/y.. = q* is obtained giving a final quotient 

q* + 2" 39 q* 

Remarks . In case the divisor is -1, the subroutine takes 

-39 
the negative of x +2 x and then jumps out. In all other cases the above 

method is followed. When y p = the division method gives an exact quotient when 

an exact quotient of 79 or fewer binary digits exists . This is particularly 

-39 
useful for integer work. The correction process works best if y +2 y is a 

standardized number » 

The mathematics of the correction process ares 

-39 -39 

Given x ± + 2 ^ x = q y + 2 •" r and y > r > 0, 

determine an integer cA such that 

x x + 2-59 X£ « ^ + * 2 -59 )(yi + 2 -59 yg) + 2 -?V - 2" 5 ^ y 2 ) 

with y > r* > . ' 
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This gives the equation r = Jl, y + q y 2 + r* . 

Since 1 > y.. > l/2 and |q | < 1 we have |q,y p | < 2y 1 and 

(-2 *A)y 1 < r - q^ - <J~ y ± < (3 - *L ) y ± or 

(-2-ol) y < r* < (3 -cO y . Since y > r* > we must have -3 < A < 3, or, 

since oC is an integer ^ = one of -2, -1, 0, 1, 2. 

Then r* = q* y^ + 2~ 59 r** with |r**| < .1/2 y . Consequently 

x x + 2~ 39 x 2 = q* (y x + 2~ 39 y 2 ) + 2~ 39 (q* y x + 2~ 59 r** + 2" 59 - r y^) 

= (qj + 2" 59 q|)(y 1 + 2" 39 y 2 ) 



+ 2"'" (- q* y 2 + p**- , y 2 ) 



■39 
Hence the maximum error in the quot ient is 7 x 2 

OVERFLOW ANALYSIS 

'1 1 A 2 



we have |-q*y 2 + r** - ^ -^ | < 6 l/2 y x < 7(y + 2* 59 y ) 

,-78 



-59 
The initial division is y into x +2 x ° It is known 



that 1 > y 1 > 1/2 while -l/2 < 1 + 2~^ 9 x < l/2 since the dividend lies in the 

floating accumulator . Consequently the first quotient is less than 1 in absolute 

value unless x = -I/2, x = and y = 1/2. In this case, the ILLIAC quotient 
_39 -L d. ± 

is -1 + 2 1 when corrected this yields q = -1 and r = 0, The correction 
, process yields nothing further since y = 0. In all other cases we have 



-39 „ ^ „ ^ iv , -39 ,_ 



y-L + 2 y 2 > y x > \\ + 2 ^ x 



Therefore q of the equation 

-59 -59 

x 1 + 2 Jy x 2 = q^ + 2 Jy r y ± > r > 

satisfies -1 < q < 1. The correction ck satisfies the relation 

r = <k y x + q.^ + r* with y > r* > . 
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Consequently if q > we must have c^< since y > 0, y > 0, r* > and 
y > r. If q < 0, then, as < r* < y , should we have cX, < - \ 3 then 
r* +ok y < and q y + <^ y + r* < 0. Contradiction as r > 0, Hence it has 
been shown that |q* I = Iq.. + o(. I < q_. I and overflow cannot result. 

SIGN CHANGE SUBROUTINE 

This little subroutine takes the negative of the double 
precision fixed point number in S3,, IS 3 and puts it back in S3, 1S3. When the 
subroutine is left, A is cleared to zero. Enter by 

P F5 pF 
p + 1 26 27% 

Remember that -(-l) = -1 in ILLIAC. 

FIXED POINT INTEGER-FRACTION READ 

Read one unsigned integer or fraction of 12 or less decimal 
digits from tape, convert the number to its binary equivalent 5 and store it at a 
specified address. 

For integer input, make the entry 

p JO nF 50 pF 
p + 1 26 327L 

The subroutine reads the decimal digits of a unsigned integer counting the number 
of digits read* Input is stopped by any character K, S, N, J, F, L or any 5 hole 
character other than a space. If c is in A on termination and c > then c-10 is 



in location when the routine is left. If c < then c-15 is in location 0. 

k-1 
If k digits are read, Q, contains 5 x 10 on exit fr 

the routine. The integer is stored at location n of the Williams Memory . 
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For unsigned fraction input make the entry 

p 50 nF 50 pF 
p + 1 26 J27L 

The subroutine reads the decimal fraction and converts it to a binary 
fraction using essentially the techniques of library routine N12 . The binary 
fraction is stored at location n. The input is terminated in the same way as 
for integer input. The contents of location are, on exit, C-10 or C-15 
according as C > or C < 0, See above. 

Caution '. This subroutine uses a 91 ^F instruction to read from tape. 
The routine is sensitive to all 5 hole characters except the space. The sub- 
routine is programmed to skip spaces . 

Other parts of routine A7 use 80 IF instructions. Indiscriminate use 
of the 80 IF instruction between two 91 ^F instructions can cause ILLIAC to read 
the same 5 hole character twice rather than advance the tape in the reader. The 
programmer who wishes to use this subroutine must be prepared to program around 
this difficulty. 

PRIM 1 ONE k-DIGIT POSITIVE INTEGER 

-59 
This subroutine should be entered with the positive integer n x 2 

in Ao There are four entries corresponding to controls of zero suppression 

and spacing «, The spacing is that specified by the last interpretive 9F nF 

instruction executed. The integer may contain k < 12 decimal digits. 

p J2 kF 50 pF Punch or print a k digit integer with- 

p + 1 26 5^9L out spacing or zero suppression on the left. 

p JO kF 50 pF Punch or print a k digit integer with a space 
n 4- 1 26 ^-QL after each n decimal digits but without zero 

' suppression on the left. 
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p 52 kF 50 pF Punch or print a k digit integer without 

^r ■,»«,. spacing but with zero suppression on the 

P + 1 26 3^9L ^ ft> 

p 50 kF 50 pF Punch or print a k digit integer with a 
ofs *koT space after each n digits and with zero 

p + ^ ■* suppression on the left. 



k 
Caution s if k digits are requested but the integer n > 10 _, this subroutine will 

not punch or print the extra digits . 

METHOD OF EXPONENT CONVERSION 

The method of exponent conversion used is based on the observation that 
the fractions l/2, 10/l6, l(r/2 > 10 / 2 > 10 / 2 can he represented exactly 
as single precision binary fractions. Since these fractions are nearly. equal to 1, the 
amount of standardization necessary for exponent conversion is small. 

The process of obtaining the floating binary equivalent of, for example, 
10 , is that of forming a double precision product 10 ^ J'P x 10 /2 x l(r/2 
and assigning the binary exponent 80. For a decimal exponent E in the range 
< E < 153 we have 

< E = a n + 3a^ + 6a r + 15a, c < 153 
— x 1 3 o 15 — 

which gives the inequalities a < 10, a^ ■+ a < 2, a < 1 and a.. < 2. Thus for 
exponents in this range the worst case is E = 1^9 for which a = 9* &g = 2, 
a, = 2o 23 single precision multiplications is the maximum number required to 
determine 

E x E f 
10 = f • 2 . 

k 
Let E .= 3K+ k with < k 2 < 2 Then k^ < 50 and (l0 5 /2 10 ) 1 > • 3 and 

p ^10 ^~\ 
(10/16) « (10 D J 2 ) > • 19. Hence during input the process of converting the 

E 
x 
decimal power 10 to a binary equivalent 

E f 
f • 2 
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can be accomplished with at most two standardizing shifts of f . Further 

-E_p -E 

l/f • 2 is a binary equivalent of 10 . The problem of exponent conversion 

during input is thus reduced to the following steps s 

(1) Read fractional part and store in floating accumulator 

(2) Read decimal exponent 

Ep E 
(5) Compute binary equivalent f • 2 of 10 X and standardize f« 

(k) Form a double precision product or quotient of the floating 
accumulator with f • 2 



The problem of binary to decimal conversion during output is essentially 
the same as the above process . However, while in binary l/2, 10/l6, . . . represent 
positive powers of 10, in the reverse process they represent negative powers of 
2. That is to say 

(10 5 /2 10 x 10 6 /2 2 ° x vP/2») x 10- 24 - 2- 80 
Hence for binary exponents in the range -512 < E < we have 

-50 b 50 - 20 * 20 - 10 b 1Q - h ^ - ^ = E x 
which gives the inequalities b < 10, b 2Q + b < 2, b < 1, h. + la < 3, 

The worst case is E = -^93 requiring 25 multiplications. The worst 
scaling problem is 

(10 3 /2 10 ) 50 (10/16 f (1/2 ) 5 > .02 o 

Thus after obtaining the decimal exponent equivalent at most one standardization 
by 10 is required . 

Conversion of a positive binary exponent by means of this process requires 
a division of the fractional part by the adjusting factor computed » Prior to this 
division the fractional part is multiplied by l/lO in double precision to eliminate 
the possibility of division overflow. 
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During input the coded conversion process will convert any decimal 
exponent to a binary exponent and place this exponent in the floating accumulator . 
Converted exponents out of the range + 512 cause no difficulty as long as no 

attempt is made to store the number from the floating accumulator. The exponent 

39 
range in the accumulator is + 2 . 

During calculation of a series of products binary exponents out of 

range may be obtained. Such exponents may be converted to decimal equivalents and 

a floating decimal answer punched «, This will cause no difficulty so long as the 

computed conversion factor is 

(l0 3 /2 10 ) 10 (10/16) k (1/2) X > .01 , 

where b^ < 2 and b < 3° Failure of this scaling rule may result in a division 
overflow during conversion of positive binary exponents but will not cause a 
failure of the routine during conversion of negative binary exponents „ 

The number in the floating accumulator is not standardized before 
exponent conversion for output . Hence if the fractional part has become less 
than l/lO, one or more zeros may appear on the left of the fractional part which 
is punched. This only means that enough arithmetic was performed without 
standardization to lose digits of accuracy . 

FIXED POINT OPERATION 

The only essential difference between the fixed point and floating point 
modes of operation is in the number storage . These differences have been described. 

Fixed point addition operates exactly as floating point addition, "When 
a number is brought to the number register 3S3, kS3, 5S3 5 location 5S3 is cleared to 
zero. The number in 2S3 gives the number of previous shifts right of the floating 
accumulator to correct for overflows «, In forming the sum this number is used to 
position the summands . Thus an arithmetically correct result may be constructed 
using the count in 2S3° 

When operating in the fixed point mode the storage subroutine does not 
take into account location 2S3° The programmer must write his program in such a 
way that if overflow occurs it is detected by, for example, testing 2S3 for zero. 
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The division subroutine requires the divisor to be standardized- In 
floating point operation the divisor is brought from the memory and hence is 
standardized. In fixed point operation this need no longer be the caseo 

Use of the standardize instruction N2 F will standardize the floating 
accumulator „ If k shifts are required then N(2S3) r = N(2S3) - ko This is 
valid for both modes of operation. 

LOCATION OF PARTS OF A? 

OL Interpretation of instructions 

17L Entry for next interpretive instruction 

19L Make the reference to a b -register 

22L Switch for interpretation of 2 function digit 

38L b -registers 

k6h Fetch one number from the memory 

230L Multiplication subroutine 

244L Division subroutine 

27% Take the negative of the floating accumulator 

284L Standardize subroutine 

327L Integer-fraction input subroutine 

3^9L Print one integer subroutine 

HIDDEN COKSTAMS 



15% 00 F 00 521F 

385L 2^/l0 n (single precision) 

kOkL 00 F 00 10F 

50% l/lO (double precision) 



DAT E July 21, 1958RTS 5/18/5 9 
PROGRAMMlID B Y *?+* , , U %** 
APPROVED B Y § ,$* ]%jjbuL^J 



lgr 



LOCATION 



ORDER 

"OOT(E7T 
00 59F 

LO 12L 

LO 13L 

ifO 2L 

00 F 

00 F 
kO 2F 
50 2F 
36 1% 
L4 392L 
36 1% 

01 7F 
I/f 7L 
te 8L 
k-2 h6L 
S3 22L 
^-0 % 
26 F 
00 F 
00 F 
kO F 
00 F 

LL 4095F 
LL 30?2F 
LL k09k-F 
it-K 4076F 
5S F 
S3 F 
5S IF 
S5 F 
50 F 
S5 20F 
00 IF 
00 2F 



NOTES 
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A 7 



7 
8 

9 
10 

11 
12 

13 
ik 

15 
16 



Bring next interpretive 

from memory . 

Closed subroutine entry 



instruction 



Jump to function switch 



Fixed point -floating point switch 



Constants 



LOCATION 


ORDER 


NOTES PAGE 


2 


IT 


36 1L 




Next interpretive Instruction 




18 


L4 i4l 

22 1L 




after 8 J 




19 


01 3F 

Lh 5IL 




Ma.Ke b - reference 




20 


42 50L 
42 181L 








21 


01 k¥ 

26 6l 








22 


26 kbh 
00 t4l 




bO Function Switch 




23 


2b 401* 
00 TIL 




bl 




2k 


Id 13L 
26 80L 




b2 




25 


27 80L 
00 F 




b3 




26 


2,6 46L 
00 101L 




b4 




- 2? 


26 H6l 
00 125'L 




b5 




28 


26 46L 
00 129L 




b6 




29 


26 k6L 
00 135L 




b7 




30 


81 4f 
26 414L 




b8 




31 


L5 21% 
26 46% 




b9 




32 


41 F 
26 14% 




bK 




33 


L5 S3 




bS 






26 160L 


.. 





A 7 



LOCATION 


ORDER 


NOTES 


PAGE 5 


54 


26 46l 

00 201L 




bN 




55 


L5 2F 
26 205L 




bj 




56 


L5 2F 
26 215L 




bF 




57 


50 9L 
26 225L 




bL 




58 


00 F 
00 F 




"b -registers 




59 


00 F 
00 F 








4o 


00 F 
00 F 








in 


00 F 
00 F 








42 


00 F 
00 F 








^5 


00 F 
00 F 








44 


00 F 
00 F 








45 


00 F 
00 F 








46 


50 9L 

L5 F 




Bring a number 


from the memory 


47 


42 6lL 
01 12F 








48 


40 F 
L5 2F 








49 


52 50L 
LO 595L 
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50 


36 51L 
50 F 






51 


JO 396L 
S5 38L 






52 


Lk F 
^2 55L 






53 


lA 68L 
*K) 62L 






5^ 


LO 6% 
36 62L 






55 


ii-1 5S3 

L5 F 






56 


ho 3S3 
F5 55L 






57 


k-2 58L 
LO 70L 






58 


32 63L 
50 F 






59 


L3 10L 
32 6lL 






6o 


S5 F 
JO 11L 






6i 


k2 5S3 








26 F 


j 


62 


83 11F 


j 




00 F 






63 


ho 3S3 
F5 62L 






6k 


4o 65L 

50 F 






65 


85 11F 
00 F 






66 


kO kSJ 
50 kS3 
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67 


kl 5S3 
26 59L 






68 


85 11F 








00 1536F 




/ 


69 


85 11F 

00 2560F 






70 


kl 5S3 
L5 1024F 






71 


F5 27L 








k2 79L 




Execute b0 5 bl inst. 


72 


L5 3S3 
10 IF 






73 


kO 3S3 
26 75L 






7^ 


F5 26L" 
22 71L 






75 


SI F 
kO 4S3 






76 


50 kS3 
32 78L 






77 


JO 59OL 
Fl 3S3 






78 


26 79L 
Ll 3S3 






79 


40 3S3 
22 F 






80 


50 2F 








iA 15L 




Execute b2, b3 insto 


81 


kO F 

01 3F 






82 


IA 83L 
k2 85L 




Set switch for 8, 9, K, S, N 


83 


10 3F 
S5 86L 


i 
i 

j 
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Qk 


h6 F 
32 96L 




Jump for b < 7 




85 


L5 S3 
26 F 




Switch 




86 


36 87L 
26 17L 




b = 8 




87 


L5 F 
22 1L 




b = 9 




88 


hi 10L 
26 17L 




b = K 




89 


49 10L 
26 17L 




b = S 




90 


L5 S3 
40 3S3 




b = N, Execute standardize N2 .inst- 
Entry after 8j to standardize FA. 




91 


50 1S3 
F5 91L 








92 


26 2QhL 
10 IF 








93 


40 S3 

S5 F. 








9h 


ho 1S3 
F5 2S5 








95 


LO 5S3 
40 2S3 








96 


26 17L 
L5 50L 




Execute count-jump on b-register 




97 


h2 98L 
hz 99L 








98 


L5 16L 
ik f 








99 


50 F 
40 F 
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100 


36 17L 
26 87L 






101 


L5 3S3 
10 IF 




Execute b4 instruction 


102 


40 3S3 
L5 5S3 






103 


LO 2S3 
40 F 




E - E 
y x 


104 


L3 F 
32 116L 




IE - E I = 0. Add immediately 
1 y x 1 J 


105 


LI F 
32 111L 






106 


t c crt-z 

40 2S3 




E > E • Interchange 
y x 


107 


S5 F 
50 1S3 






108 


40 1S3 
L5 S3 






109 


40 4S3 
L5 3S3 






110 


4o S3 
L5 4S3 






111 


4o 3S3 




Set for positioning. 




L7 F 




- • ■ 


112 


^2 115L 
LO 393L 




|E - E I - 79 > Skip addition 
1 y x 1 — * 


113 


36 17L 
L4 394L 




|E - E i - 64 > Too many shifts 
1 y x 1 - J 


114 


32 122L 
26 115L 






115 


L5 3S3 
10 F 




Position and add. 


116 


40 3S3 
L5 1S3 
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117 


Sk F 
kO 1S3 








118 


50 1S3 
32 120L 








119 


JO 390L 
F5 3S3 








120 


26 121L 
L5 3S3 








121 


Lk S3 
kO S3 








122 


22 139L 
Lk 389L 




In case > 6k shifts 




123 


k2 115L 
L5 3S3 








22k 


10 63F 










22 11% 








125 


L5 3S3 
10 IF 




Execute b5 instruction 




126 


kO S3 
^0 S3 . 








127 


L5 5S3 
kO 2S3 








128 


22 139L 
00 F 








' 129 


L3 10L 
32 130L 




Execute "b6 instruction 




130 


L5 387L 
IA 2S3 








131 


LO 5S3 
kO 2S3 








132 


22 132L 
F5 132L 








133 


26 2^L 
50 1S3 
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l?k 


22 139L 
00 F 






135 


L3 10L 
32 I36L 




Execute JoJ instruction 


136 


LI 387L 
iA 2S3 






137 


L4 5S3 
kO 2S3 






138 


L5 ^S3 
F5 138L, 






139 


26 230L 
LJ S3 




Entry to correct FA for overflow 


i4o 


36 143L 
F5 2S3 






l4l 


kO 2S3 

L5 S3 






lk2 


10 IF 
kO S3 






lA3 


S5 F 
kO 1S3 






Ikk 


26 17L 
00 F 






1^5 


L5 2F 
46 F 




Execute bK instruction 


lk6 


32 155L 








L3 F 


- 




ikj 


36 15^L 
L5 F 






148 


50 386L 
00 10F 






1^9 


kO 3S3 
F5 149L 







A 7 



LOCATION 



ORDER 



NOTES 



PAGE 10 



A 7 



150 
151 
152 
153 

155 
156 

157 
158 

159 
160 
161 
162 
I63 
164 
165 
166 



26 284l 
10 IF 
kO S3 
hi 1S3 
L5 159L 
LO 5S3 
kO 2S3 
26 17L 
kl S3 
kl 1S3 
26 153L 
L5 50L 
42 157L 
50 386L 
LI F 
40 F" 
26 17L 
00 F 
00 F 
00 521F 
kO 3S3 
50 1S3 
L3 10L 
36 199L 
50 1S3 
F5 162L 
26 284L 
40 3S3 
L3 5S3 
36 174L 
L5 387L 
Sk F 
40 4S3 
50 4S3 



Entry to clear floating accumulator 



Execute bS instruction 



Jump for fixed point 



Standardize 



test for zero. 



Round 
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I67 

168 
I69 
170 

171 
172 

173 
17^ 
175 
176 

177 

178 

179 
180 
181 
182 
I83 



32 172L 
JO 390L 
F5 3S3 
ho 3S3 
32 172L 
LI 3S3 
32 172L 
^9 3S3 
F5 389L 
L4 2S3 
26 173L 
F5 2S3 
LO 5S3 
36 175L 
4l 3S3 
50 386L 
JO 11L 
Sk F 

ho ksj 
SO F 
LO 398L 
32 195L 
30 9L 
01 12F 
kO F 
L5 2F 
32 181L 
LO 395L 
36 182L 
50 F 

JO 396L 

S5 F 
Lk F 
k2 186L 



Correct exponent for stand 



Set no . = » 



Pack exponent and test for overflow 



Unpack address 



Determine b-modif ication 



B-modify address and set store 
instructions , 



LOCATION 



184 

185 

186 

I87 

188 

I89 

190 

191 

192 

193 

19^ 

195 

196 

197 
198 

199 



ORDER 



Lk I96L 
kO 191L 
LO 197L 
32 190L 

L5 3S3 
Aj-0 F 

F5 186l 
42 189L 
LO I98L 
36 192L 
L5 4S3 
h-0 F 
26 I7L 
L5 3S3 
86 11F 
00 F 
F5 191L 
4o 194l 

50 F 
L5 4S3 
86 11F 
00 F 
26 17L 
FF 57F 
86 iif 
00 1536F 

86 11F 
00 256OF 
L5 3S3 
ho 102 kF 

L5 S3 
00 IF 
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W.M. or drum? 



FF stop for overflow, 
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200 


40 5S5 
25 17% 






201 


L5 S5 
52 2.05L 




Execute bN 


202 


22 202L 
F5 202L 






205 


26 279L 
L5. 5S5 






20k 


56 7^L 
26 101L 






205 


52 206L 
k6 206L 


, 


1 

1 


206 


26 F 




Execute bj 




k6 207L 






207 


LI F 
kO F 




i 
1 


208 


L5 50L 
h-2 210L 


i 


1 

i 

i 


209 


k2 215L 








50 596L 


i 


! 


210 


JO F 


I 


1 




L5 F 


\ 


I 


211 


Sk F 
ho F 


1 


; 


212 


50 F 
JO 599L 




! 


215 


S5 F 
40 F 






214 


26 17L 
00 F 






215 


52 219L 
L^ 592L 




Execute bF 


216 


52 221L 








h6 218l 


1 
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217 
218 
219 
220 
221 
222 



224 
225 
226 

227 
228 

229 
250 
231 
232 
233 



Jj-7 219L 
92 131F 
50 F 
26 17L 
00 F 
46 220L 
L5 F 
22 207L 

00 F 

46 211L 

47 221L 
26 17L 

01 12F 
40 F 

47 207L 
L5 2F 
36 208L 
L4 392L 
36 207L 
L5 50L 
42 228L 

42 213L 
50 397L 
JO F 

43 210L 
22 210L 
42 238L 
L4 397L 
46 234l 
46 235L 
LJ 388L 
74 S3 
40 IF 
S5 F 



Execute bL 



Multiply subroutine 
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254 


50 F 
74 1S5 




■ 


255 


50 F 
74 S5 






256 


40 S5 








L5 IF 






257 


S4 F 
4o 1S5 






258 


50 1S5 








52 F 




Link 


259 


JO 590L 
L5 IF 






24o 


56 242L 
Fl 586L 






24l 


L4 S5 








22 242L 






242 


F5 S3 








40 S5 






243 


27 237L 
00 F 






244 


42 276L 








L5 3S5 




DiTide subroutine 


2^5 


40 F 
32 248L 




Save Sign 


246 


L5 255L 
42 79L 






247 


26 75L 




Absolute value of divisor 




52 248L 




l y l + 2 " 39 y 2 1 


248 


26 277L 
S5 F 




Jump for -1 


249 


40 4S5 
50 1S5 
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250 


jo 58% 




Save do 






S5 276L 






251 


kO IF 
50 1S5 








252 


L5 S5 
66 5S5 




S"fc 

1 division 




255 


kO 1S5 
S5 2kjL 




Store remainder r. 




25^ 


kO S5 
56 2.66L 




Store quotient q.. 




255 


Fl 586L 
Lk IF 








256 


Lk 5S5 
lA 1S5 








257 


kO 1S5 
52 260L 




Obtain correct remainder 




258 


Lk 5S5 
40 1S5 








259 


L5 S5 
FO 586L 








260 


^0 S5 
50 te3 




- 1 X x y 2 




26l 


71 S5 
52 267L 








262 


Lk 1S5 
52 271L 








265 


Lk 5S5 
kO 1S5 




Adjust first remainder 
for double length divisor 




2.6k 


L5 S5 
FO . 586L 








265 


kO S5 
27 262L 








266 


L5 IF 
LO 5S5 
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267 


22 256L 
Lk 1S5 








268 


52 271L 
kO 1S3 








269 


F5 S5 
kO S5 








270 


L5 1S5 
LO 5S5 








271 


26 268L 
LO 5S5 








272 


56 278L 
80 IF 








275 


L4 5S5 
50 5S5 




2 division 




2Jk 


66 3S5 
10 IF 








275 


SJ F 
kO 1S5 








276 


L5 F 
52 F 




Link 




277 


L5 25OL 
26 279L 




Change sign of quotient when 
divisor negative 


! 


278 


L^ 5S5 
22 268L 






! 


279 


k2 285L 
LI 1S5 




Take negative of FA 




280 


52 282L 
FO 590L 




Subroutine 


■ ' 


281 


kO 1S5 
Fl S5 








282 


26 283L 
LI S5 








285 


ho S5 
25 F 




Link 
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284 


h2 322L 
kl 5S3 




Standardize subroutine 




285 


L5 3S3 
32 286L 








286 


26 29IL 
SI 8F 




Form - |x + 2 ~ 59 x | 




287 


36 290L 
kO H-S3 








288 


50 ^S3 
JO 390L 








289 


Fl 3S3 
26 291L 








290 


LI 3S3 
32 322L 








291 


00 IF 
36 303L 




1 
"Inner loop" 




292 


00 IF 
36 304L 




Shift until sign of A changes. 




293 


00 IF 
32 305L 








29^ 


00 IF 
36 307L 








295 


00 IF 
32 308L 








296 


00 IF 
32 309L 








297 


00 IF 
36 3HL 








298 


00 IF 
32 312L 








299 


00 IF 
36 31^L 
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40 4S3 




L5 5S3 


301 


iA 286L 




k2. 5S3 


302 


L5 hS3 




26 292L 


303 


40 kS3 




22 316L 


30^ 


ho kS3 




L5 389L 


305 


22 315L 




hO kSJ 


306 


F5 389L 




22 315L 


307 


ho hs3 




L5 391L 


308 


22 315L 




40 4S3 


309 


27 315L 




40 4S3 


310 


L5 389L 




26 315L 


311 


40 4S3 




F5 389L 


312 


26 315L 




40 4S3 


313 


L5 391L 




26 315L 


31^ 


40 4S3 




F5 391L 


315 


F^ 391L 




Lh 5S3 


316 


^2 5S3 




L5 3S3 
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Count 8 and repeat above 



Count 0, 1, . .., 7 



According to the jump executed 



if x > 0, take negative again 
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317 


32 318L 
L5 ^S3 






318 


22 321L 
SI F 






319 


36 323L 
kO 3S3 






320 


50 3S3 
JO 390L 






321 


Fl 4S3 
10 IF 






322 


FO 390L 
22 F 




Link 




LI kS3 
32 32^L 






324 


22 321L 
Fl 386L 






325 


Lk 5S3 
k2 5S3 






326 


2S 322L 
00 F 






327 

1 


K5 F 
k2 3^8l 




Integer-Fraction input subroutine 


328 


k6 3^8L 
36 330L 






329 


L5 327L 
22 330L 






330 


F5 330L 
42 3^5L 






331 


41 kS3 
50 405L 






i 532 


26 337L 
LO 4o4L 




Read loop 


333 


32 338L 
10 3F 
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33^ 
335 
336 

337 
338 
339 
3^0 
3^1 
3^2 

3^3 
3^4 

3^5 
3^6 

3^7 
3^8 

3^9 



F4 4S3 
00 2F 
Fk 4S3 
00 IF 
40 4S3 
00 IF 
91 4f 
32 332L 
LO 39^L 
ho F 
L7 F 
36 3^1L 

00 4f 
26 337L 
S5 F 

10 h-F 

01 kF 

42 3kkL 
k2 346L 

L5 ^S3 
50 F 
22 3^5L 
26 F 
SO F 
66 F 
10 IF 
SJ 400L 
kO F 
22 F 
40 F 
kl 2F 



Test for spaces 



Decode for 5 x 10 



Divide for fraction input 



Link - Store Number 
Integer Print Routine 
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350 


K5 37^L 
42 375L 




Plant link. 




351 


46 2F 
32 352L 




Store no. of digits 




352 


27 353L 
L5 365L 








353 


42 353L 
S5 F 




Set zero suppression 




354 


00 6f 
32 356L 








355 


L5 350L 
46 370L 




Set spacing 




356 


22 357L 
L5 365L 








357 


46 370L 
L5 384l 




Set count on number of digits 




358 


LO 2F 
40 2F 








359 


L5 f 
50 385L 








360 


32 361L 

LI 364L 




Put rounding constant in A. 




361 


26 362L 
L5 364l 








362 


74 F 
36 364l 








363 


l4 385L 
L4 385L 




x 2 55 /l0 1] - 




364 


10 35F 

40 F 




Store one character at 




365 


S5 377L 

40 IF 




Store residue at 1 




366 


L5 2F 
36 382L 
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567 


L5 F 
52 570L 




Test for zero . 


368 


h3 555L 
L5 F 






569 


00 56F 

82 k-F 




Print one character. 


570 


22 F 
L5 555L 




Zero suppress? 


571 


10 IF 
SJ F 






572 


52 568L 
F5 2F 






575 


00 20F 
52 568L 




Last character to be printed? 


57^ 


26 58IL 
F5 2F 




Count 


575 


kO 2F 
52 F 




Link. 


576 


50 IF 
75 kokh 






577 


22 56^ 
L5 221L 




Space between groups of digits? 


578 


lA 597L 

k6 221L 






579 


52 57^L 
L5 211L 






580 


L4 597L 
46 221L 






581 


92 965F 
22 57^L 




Print one space. 


582 


Fk 597L 
kO 2F 






585 


26 57^L 
00 F 
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384 


00 HP 
LL 4084F 




Count. 


385 


2S 4015F 
LN 755F 




235/xo 11 


386 


00 F 








00 F 




Constants 


387 


00 F 
00 512F 






388 


20 F 
00 F 






389 


00 F 
00 IF 






390 


7L 4095F 
LL 4095F 






391 


80 F 
00 3F 






392 


70 F 
00 F 






393 


00 F 
00 79F 






39^ 


00 F 
00 15F 






395 


10 F 
00 F 






396 


00 F 

00 4095F 






397 


LL 4095F 
00 F 






398 


00 F 

00 102^F 






399 


LL 4095F 
00 4095F 
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400 


74 1701F 
28 2048F 




"11 k -^Q 
5 x 10 ^ Table of 5 x 10 x 2 Jy 


401 


OS 262 7F 
S7 1024F 




5 x 10 


402 


00 47F 
KL 128F 




5 x lo" 7 


403 


00 4f 
n4 288OF 




5 x 10 


4o4 


00 F 
00 10F 




-39 
10 x 2 Jy 


405 


23 619F 

j6 1616F 




Key word 


4o6 


00 476F 
j6 128OF 




5 x 10 


407 


00 F 
ON 848F 




4 
5 x 10 


4o8 


01 672F 

5L 512F 




5 x 10 9 


409 


00 F 
00 50F 




5 x 10 


4io 


00 F 
7K 288F 




5 x 10 5 


4li 


00 F 
00 500F 




2 
5 x 10 


412 


00 F 




5 




00 5F 






413 


00 F 

01 904f 




5 x 10 5 


414 


lo 4o4l 




Execute 88 instruction 




40 5S3 




Store sign 


415 


36 4i6l 








FF 58F 




Format error. 


4i6 


JO S3 
50 4i6l 




Read first part of number as integer 
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4l7 


26 527L 
S5 551F 






4l8 


40 5S5 
L5 F 




Store 5 x 10 x 2 Jy 


4-19 


L4 594l 
52 458L 




Test termination for "." 


420 


lo. 4o4l 

40 F 






421 


L7 F 
56 440L 




Jump for floating decimal format. 


422 


50 1S5 








50 422L 




Read fractional part 


425 


du Jd. { J_i 

L5 4l7L 




insert correct exponent. 


424 


42 2S5 
50 1S5 






425 


L5 S5 
52 429L 






426 


10 IF 




Test integer part for overflow 




FO 590L 




and correct. 


427 


40 S5 
F5 2S5 






428 


40 2S5 
S5 595L 






429 


40 1S5 
L5 5S5 




Take negative of ab» value if minus. 


450 


52 451L 
F5 450L 






451 


26 279L 
L5 S5 






452 


40 5S5 
50 597L 


' 




455 


50 1S5 
F5 455L 




Standardize 
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434 


26 284L 
10 IF 






435 


40 S3 
S5 F 






436 


40 1S3 
L5 2S3 






437 


LO 5S3 
42 2S3 






438 


26 17L 
4l 1S3 




Set 1S3 = when < 12 digits 


^39 


22 441L 
00 F 






44o 


50 1S3 
50 440L 




Read floating decimal numbers . 


hhi 


26 32 7L 
L5 F 




Read 2 part of number as a fraction 


442 


36 443L 
22 4l5L 




Format error if number in 3 parts « 


~~ 443 


42 455L 

4l 4S3 




Store sign of exponent. 


444 


50 1S3 
L5 S3 






445 


10 IF 

32 446l 






446 


FO 390L 




Scale floating accumulator to correct 




10 IF 




range » 


447 


40 S3 

S5 F 






448 


40 1S3 
50 F 






449 


50 4S3 
F5 449L 




k -39 
Double length division by 5 x 10 x 2 
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450 


26 244L 




Take negative for negative numbers 




L3 5S3 




■ 


451 


' 32 452L 








F5 451L 






452 


1 26 279L 








50 F 




. 


453 


JO F 




Read exponent. 


: 

! 


50 453L 






454 


26 327L 




: 


" 


LI F 






455 


40 IF 


■ 


Obtain signed exponent. 


• 


L5 F 


■ 




: 456 


40 F 








L3 F i 

1 


! 


' 


^57 


j 36 467L 






; 


F5 457L 






458 


26 564l 




Entry to exponent conversion subroutine 




L5 5S3 






^59 


40 IF 








F5 459L 






460 


26 284l 




1 
Standardize conversion factor 




40 3S3 






461 


L5 F 








32 464l 






462 


Ll IF 
L4 5S3 






463 


L4 387L 




Multiply or divide according to 




40 2S3 




sign of exponent* 


464 


22 132L 
L5 IF 






465 


LO 5S3 
l4 387L 






466 


40 2S3 
26 I38L 
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467 


L5 587L 
4b 2S3 






468 


50 1S5 
22 139L 






469 


L4 397L 
46 219L 




LF-CR? Execute 89 instruction. 


470 


56 475L 
92 151F 






471 


92 515F 
L5 218L 




Print LF-CR and delay. 


kj2 


L4 397L 




Reset spacing control. 




46 219L 






475 


L5 211L 
L4 397L 






474 


46 221L 








L5 S5 






475 


56 477L 
92 706F 




Print + sign 


476 


L5 519L 
26 279L 




Take absolute value 


477 


92 642F 

92 963F 






478 


L5 587L 
LO 2S5 






479 


40 F 
32 481L 




512 x 2~ 59 - (E + 512) o2' 59 

.A. 


480 


50 506L 




Multiply fractional part "by l/lO 




F5 480L 




for later division 


481 


26 250L 
L5 F 






482 


56 507L 
F5 482L 






485 


26 566L 




Compute exponent conversion factor. 




L5 5S5 




, 
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484 


40 F 








LO 505L 






485 


36 488L 








F5 48% 






486 


: 26 586L 


■ 


Multiply conversion factor "by 10. 


j 


, ^0 3S3 






m 


F5 5S3 


\ 




5 


40 5S3 


i 




488 


L5 387L 








LO 2S3 






489 


36 502L 


.• 


Jump when E < 


- 


F5 5S3 






490 


, 40 2S3 


i 


Divide "by correction factor. 


■ 


F5 490L 


f 




491 


26 244l 




j 


; 


50 1S3 






492. 


L5 S3 


| 


i 




00 IF 




| 


493 


40 F 








4o S3 \ 




i 

1 


494 


S5 F 




; 




40 1S3 ! 




i 


495 


L5 F ' 




Multiply "by 10 if possible * 




F0 505L i 




i 


496 j 


36 509L j 




i 


1 


F5 496L 


j 


j 

3 


497 \ 


26 586L \ 
32 498L 




Check if xlO produces overflow. 


498 


26 509L 








40 S3 






499 


S5 F 






nMwun-, ■ , -im-m^mv 


40 1S3 
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500 


L5 2S3 
FO 386L 






501 


kO 2S3 
26 509L 






502 


LI 5S3 




Multiply by conversion factor if 


: 


kO 2S3 




E < 0. 

X 


503 


: L5 ^S3 








F5 503L 




: 


50^ 


2.6 230L 








22 ^91L 






505 


ON 32T6F 








NN 3276F 




l/lO double precision 


506 


66 1638F 
66 1638F 






507 


kO 2S3 
LJ S3 




1 


508 


32 491L 
22 5i«-7L 






509 


kl F 








L5 2F 






510 

1 


k6 F 




Determine number k of digits 




L5 F 




to print 


511 


LO 384L 




jump if > 12 digits* 




36 531L 






512 


L5 F 
k6 51**-L 






513 


\6 52OL 
50 F 






5lA 


50 F 
50 51^L 






515 


26 551L 
LJ 183 




Multiply by 10 , rounded * 


516 


32 517L 




Round product to nearest integer 




F5 S3 


' Test rounded product for overflow. 
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517 


kO S5 
L5 S5 






518 


LO F 
52 5^7L 






519 


L5 S5 








50 ^77L 






520 


JO F 
50 520L 




Print k digit integer 


521 


26 5^9L 








92 965F 




Space 


522 


L5 2S5 
56 525L 






525 


92 706F 
Ll 2S5 




Print sign of exponent. 


524 


40 2S5 
22 525L 






525 


92 6H-2F 
L5 2S5 




Print exponent » 


526 


52 5F 
50 526L 






527 


26 5^9L 








L5 219L 






528 


L^ 597L 
52 529L 






529 


26 17L 
92 IF 






550 


92 967F 




Exit. Print delay and 2 spaces if 




26 17L 




LF-CR will not preceed next number. 


551 


LO 597L 
k6 552L 




when > 12 digits, take 

in k-ll 
x 10 


552 


50 F 
50 552L 






555 


26 55IL 
50 ^OIL 


■ 
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534 


00 IF 
7J 1S3 




and then x 10 for residue 


535 


40 1S3 








LO 401L 




'■ 


536 


LO 401L 
36 544L 




Test low order part for overflow. 


537 


22 5^5L 
L5 532L 






538 


46 539L 
L5 S3 




Print higher part 


539 


JO F 
50 539L 






540 


26 3^9L 
92 IF 




delay 


54l 


L5 1S3 
50 F 


' 


. ! 


542 


JO 11F 
50 542L 




Print lower part. 


5^3 


26 3^9L 
22 521L 






544 


41 1S3 




Correct high order part and 




F5 S3 




test for overflow. 


5^5 


40 S3 
L5 S3 






5^6 


LO F 
32 5^7L 






5^7 


22 537L 








F5 2S3 




When rounded product overflows, 


548 


40 2S3 




replace "by l/lO double precision 




L5 505L 




and repeat. 


549 


40 S3 
L5 506l 
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550 


40 1S3 






551 


26 509L 
4l F 


! 


k-1 -59 
Given k, determine 5 x 10 x 2 


s 


K5 F 




■ 

1 


552 


; k2 563L 




■ 


\ 553 


k6 F 
L5 F 
| 00 IF 




-59 

from table multiply x.. ' + 2 x 


' 554 


! L4 f 


; 




f 


L4 558L 




i 

' ' j 


555 


; 46 556L 


■ ; 


| byio k 




L5 405L 






556 


: 10 F 








01 kF 






557 


Lk 3kjL 




1 • 




42 558L 






558 


50 8F 
L5 F 


5 




559 


80. IF 






• 


kO F 






560 


; 50 1S3 

7 J F 






561 


, 50 S3 








74 F 






562 


kO S3 
S5 F 






563 


40 1S3 








22 F 






564 


42 58IL 




Convert Exponent subroutine 


! 


* F5 4a8L 






565 


42 573L 
22 567L 




. 


566 


42 58IL 






' i-t ■ •■ vm ■"" » - 1 


L5 428L 




■ -■ ■ . -V w,i ■ ■ -'"-' 
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575 
576 
577 
578 
579 
580 
581 
582 
583 
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^2 573L 

^9 3S3 
hi kS3 

L7 F 
50 386L 
00 20F 
kO 5S3 
L3 E 
32 58IL 
L5 ^32L 
k2 583L 
k2 584L 
50 ^S3 
L5 F 
iA 5S3 
36 583L 
L5 583L 
1,4 391L 
k2 583L 
42 58^L 
L5 573L 
Lk 391L 
^2 573L 
LO 582L 
32 573L 
L5 3S3 
00 IF 
ho 3S3 
22 58IL 
22 F 
50 hS3 
L5 6lOL 
40 5S3 
7 J F 
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584 


50 3S3 








74 F 






585 


40 3S3 








22 573L 




■■ 


586 


42 594L 








S5 F 




Multiply by 10 subroutine 


587 


^0 IF 
L5 F 






588 


00 2F 




" 




L4 F 






589 


40 F 








S5 F 




•: 


590 


L4 IF 








40 IF 




■■' 


591 


50 IF 
32 593L 






592 


JO 390L 
F5 F 


■ 




593 


26 594L 
L5 F 






59^ 


00 IF 
22 F 


' ; 




595 


LL 4046F 








00 15F 




Table for exponent conversion 


596 


LL 4081F 
00 50F 




\ 


597 


71 2813F 




10 15 /2 50 




49 2255F 


1 


598 


LL 4076F 


I 




00 6f 




■ 


599 


LL 4090F 








00 20F 
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600 


7K 288F 
00 F 




10 6 /2 2 ° 


601 


LL ko86F 
00 3F 






602 


LL ^093F 
00 10F 






603 

I 


7J F 
00 F 




10 5 /2 10 


60k 


LL 4092F 
00 IF 






605 


LL 4095F 
00 kF 






606 


50 F 
00 F 


. 


10/2 4 


607 


LL 4095F 
00 F 






608 


LL 4095F 
00 IF 






609 


kO F 
00 F 




1/2 










i 
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